Magnetic and Electric Phase Control in Epitaxial EuTiOs from First Principles 
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We propose a design strategy - based on the coupling of spins, optical phonons, and strain - for 
systems in which magnetic (electric) phase control can be achieved by an applied electric (magnetic) 
field. Using first-principles density-functional theory calculations, we present a realization of this 
strategy for the magnetic perovskite EuTiOs. 
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There is great interest in multiferroic materials in 
which ferroelectric (FE) and ferromagnetic (FM) order- 
ing not only coexist, but in which the electrical polariza- 
tion P and the magnetization M are large and strongly 
coupled P, 01 • One challenge in identifying strongly cou- 
pled FM-FE's that has received considerable attention 
in the past Q is the scarcity of such materials in nature, 
as most insulators (a requirement for ferroelectricity) are 
paraelectric (FE) and antiferromagnetic (AFM). With 
the recent advances in first-principles density-functional 
methods for predicting jJi] and in novel synthetic tech- 
niques for growing new FM-FE multiferroics, the focus 
has now turned to how to produce a strongcoupling be- 
tween the two distinct order parameters Recently, 
attention has focused on a mechanism in which magnetic 
order itself breaks inversion symmetry 0, Q, IS 0| ■ Based 
on this, remarkable control of the FE state by an applied 
magnetic field has been demonstrated in some rare-earth 
manganites [1(1. however, the natural scale of the 
spontaneous polarization thus induced is very small, of 
the order of nC/cm^. Furthermore, the magnetic state 
appears to be rather insensitive to an applied electric field 
for this class of materials. 

As a result, it is clearly advantageous to explore other 
possible mechanisms for strongly coupled multiferroism. 
A fruitful starting point for identifying such mechanisms 
is the observation, recently discussed by Tokura 5], that 
the basic physics of a strong M-P coupling involves a 
competition between different ordered states, e.g. be- 
tween a FM-FE state and an AFM-PE state. In this Let- 
ter, we present a new approach for designing a strongly 
coupled multiferroic in which the interplay of spins, op- 
tical phonons, and strain leads to such a competition. 

The criteria that a system must satisfy for this pro- 
posed mechanism to be realized are as follows: (1) It 
must be an AFM-PE insulator in which at least one 
infrared-active (ir) phonon is coupled to the magnetic 
order, (2) the spins in the AFM ground state should 
align with the application of a magnetic field of mod- 
est strength, (3) this alignment should decrease the fre- 
quency of the spin-coupled ir-active phonon, and, (4) the 
key to our approach, the ir-active mode of interest must 
be strongly coupled to strain. Epitaxial strain can have 
profound effects on the properties of thin films In 
our design strategy we use epitaxial strain to dial into 



the region of the phase diagram where a spin-phonon- 
driven destabilization of the lattice actually occurs. The 
FM-FE phase thus produced is a low-lying state com- 
peting with the AFM-PE ground state. As a direct re- 
sult of this competition, magnetic and electric phase con- 
trol can be achieved by an applied electric and magnetic 
field respectively. As the polarization of the low-lying 
FM-FE phase results from the freezing-in of a soft polar 
phonon triggered by the spin-phonon coupling, it is of the 
same order of magnitude as prototypical soft-mode EEs 
(/iC/cm^). In addition, as the phase-control region is 
approached from the low strain side, the magnetocapac- 
itance diverges, while on the high strain side there may 
be, in some cases, a phase boundary between the AFM- 
PE phase control region and a true equilibrium FM-FE 
phase. The use of epitaxial strain to exploit the relatively 
modest spin-phonon effects displayed by many bulk ma- 
terials provides an exciting new design strategy in the 
pursuit of strong M-P coupled multiferroics. 

In bulk, europium titanate EuTiOa is an AFM-PE 
that crystallizes in the cubic perovskite structure (space 
group Pm3m) with room temperature lattice constant 
equal to that of SrTiOa, a=3.905A. The Eu^+ moments 
(J=S'=7/2) order at TN=b.5K 

while neutron 

diffraction studies indicate G-type AFM order J^] . Both 
diffraction and local structural probes [3 (XANES) sup- 
port EuTiOs remaining cubic at all temperatures. Re- 
cently, Katsufuji and Takagi djj (KT) showed that at 
the onset of AFM order, the static dielectric constant 
eo undergoes a sharp reduction of about 10%, indicating 
a hardening of the lattice. Further, KT demonstrated 
that in an increasing magnetic field eo increases, saturat- 
ing at a field large enough to fully align the spins (ap- 
proximately 1.5 T). KT argue and, as we will see, our 
first-principles calculations confirm, this behavior is due 
to a coupling between the spins and an ir-active phonon 
of the type uj = lopm + A(Si ■ Sj) [H [li|. These bulk 
measurements show that EuTiOa satisfies the first three 
criteria, and a close structural analogy to SrTiOa sug- 
gests that it satisfies the key fourth criterion. In the 
remainder of the Letter we present a first-principles cal- 
culation of the ground-state epitaxial phase diagram for 
EuTiOa demonstrating that this system is indeed a real- 
ization of our proposed mechanism for strongly coupled 
multifcrroicity. 
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Within density-functional theory, the failure of the 
generalized gradient approximation (GGA) properly to 
capture the physics of strongly correlated systems is well 
established. A widely accepted approach beyond GGA 
is the GGA plus Hubbard U (GGA+U) method 
We perform first-principles density-functional calcula- 
tions using projector augmented-wave potentials within 
spin-polarized GGA-I-U approximation as implemented 
in the Vienna ah initio Simulation Package |2rl l2ll | with 
a plane wave cutoff of 500 eV and a 6x6x6 F-centered 
fc-point mesh. The PAW potential for Eu treated the 4f^ 
5s^ 5p^ 6s^ as valence states. All calculations were per- 
formed with coUinear spins and without LS-coupling. As 
expected for Eu^+ which lacks orbital degrees of freedom, 
inclusion of LS-coupling does not change the results. Val- 
ues of the Eu on-site Coulomb, U = 6 eV, and exchange, 
J_f/=1.0 eV, parameters were used that give a reasonable 
account of the magnetic exchange constants which were 
extracted by mapping GGA-t-U calculations of the to- 
tal energy for different spin configurations at T=0 onto a 
classical Heisenberg model; Espin—— Xlij Jij^i '^j (note, 
in our notation the energy per spin bond is 2J). Phonon 
frequencies and eigendisplacements were calculated us- 
ing the direct method where each ion was moved by 
approximately O.OlA. Born effective charge tensors were 
calculated by finite differences of the polarization using 
the modern theory of polarization |22l | as implemented 
in VASP. Before we proceed we note that calculations 
with GGAH-U overestimate lattice constant a by 
Therfore, we introduce a shift of the zero for (T33, the 
out-of-plane component of the stress, so that 1733=0 for 
the cubic structure at the experimental lattice constant, 
o-exv= 3.9A. Thus, the correct cubic structure is obtained 
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First, to examine the validity of a spin-phonon cou- 
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FIG. 1: Soft ir-active phonon frequency squared, lo'^, (cm~^) 
of paraelectric EuTiOa as a function of compressive epitaxial 
strain 77. Space group for r]^0 is P4/mmm; the phonon is 
polarized along the tetragonal axis, with symmetry label A2u- 



pling mechanism in bulk EuTi03 we performed first- 
principles density-functional calculations of the ir-active 
phonons and the Born effective charges, allowing us to 
evaluate eq of a single-domain, stoichiometric, defect-free 
crystal. To elucidate the role of magnetic order on the 
phonons and subsequently the dielectric response, we cal- 
culate the quantities just mentioned with the Eu spins in 
two different magnetic configurations: FM and G-type 
AFM. All calculations for cubic EuTi03 were carried out 
at the experimental lattice constant ao = 3.90 A. The cal- 
culated magnetic exchange constants were found to be in 
good agreement with experiment (shown in parentheses): 
Ji = -0.013 K (-0.014 K), J2 = +0.065 K (-^0.037 K); 
the positive sign for J2 is consistent with the unusual 
positive sign of the Curie- Weiss temperature for AFM 
EuTi03. The ground state is thus correctly found to be 
G-type AFM, with the FM state higher in energy by 0.2 
meV per formula unit. For the static dielectric constant, 
we find eo = 350 for EuTi03 with FM order and eg = 
280 with AFM order. These values compare well with 
the experimental values of 420 and 380 in fields of 1.5 
T and T, respectively The spin-dependence of 

the dielectric response is due entirely to the lowest-lying 
ir-activc phonon TOl, which we calculate to be at 70 
cm^^ and 77 cm^^ for FM and AFM respectively. This 
good agreement of the relative magnitudes of the val- 
ues for eo with the size of the experimentally measured 
field-induced change confirms KT's suggestion that the 
magnetocapacitive effect in bulk EuTi03 is indeed due 
to spin-phonon coupling, with the sign of the coupling A 
such that FM spin alignment reduces the stability of the 
low-frequency polar phonon. 

Next, we consider the effects of epitaxial strain. In 
perovskite titanates such as SrTi03 and BaTi03, it is 
well known both experimentally and theoretically that 
epitaxial strain couples strongly to the low lying TOl 
polar mode. For example, in SrTi03, a modest epitax- 
ial strain of less than 1% transforms the PE bulk mate- 
rial into a room-temperature FE 23], with P of nearly 
20 fiC cm~^ |2^. As the calculated eigendisplacemcnt 
pattern of this polar mode for EuTi03 is virtually iden- 
tical to that of the soft polar mode previously calcu- 
lated for SrTiOs and BaTi03, it is natural to expect 
similar strain coupling, with the spin-phonon coupling in 
EuTi03 adding richness to the phase diagram. 

We isolate the effects of epitaxial strain, as in p revious 
phenomenological and first-principles studies j2,4l26ll27l |. 
by imposing the epitaxial constraint on the lattice param- 
eters of the infinite crystal with periodic boundary con- 
ditions, corresponding to zero macroscopic electric field. 
For simplicity, in the present study we consider only com- 
pressive strain. We focus attention on the lowest-lying ir- 
active phonon TOl which as mentioned above dominates 
the spin-dependent dielectric response in cubic EuTiOs 
and any possible FE instability. In the cubic phase, TOl 
is 3-fold degenerate, where each phonon mode is polar- 
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ized along one of the three Cartesian coordinates. Under 
a biaxial compressive strain in the a-b plane, the lattice 
expands along the c-axis, conserving volume to a first ap- 
proximation. This lowers the symmetry of EuTiOs from 
cubic to tetragonal, splitting the 3-fold degenerate TOl 
mode into a two-fold E„ and one-fold A2„ degenerate 
mode polarized perpendicular and parallel to the c-axis 
respectively. Under compressive strain, the E„ modes 
harden and are not relevant to our study. 

In Fig. n we show the evolution of the lowest-lying 
ir-active phonon frequency with biaxial compressive 
strain 77 for FM and AFM EuTiOs. At 77=0, we have bulk 
EuTiOa where uja2u—'^toi = 70 cm~^ and 77 cm~^ for 
FM and AFM respectively, as we previously discussed. 
With increasing compressive strain, the phonon fre- 
quency decreases for both types of magnetic order. At 
approximately rjci—0.92%, the FM system develops a 
polar instability as evidenced by the vanishing ir-active 
phonon frequency. On the other hand, the AFM system 
remains stable up to approximately 77c2=l-2%. Thus, 
in the intermediate strain regime between 77R:!0.92% and 
ry«1.2%, EuTiOs has a low-lying FM-FE state in compe- 
tition with a AFM-PE state, exactly the conditions which 
have been identified as facilitating strong M-P coupling. 

It is this difference in the EE instability for the AFM 
and FM orderings that leads to a giant magnetocapac- 
itive effect near the critical strains, 77c- To see this, in 
Fig. Ela) we show the component of the static dielectric 
tensor along the c-axis, ess, calculated from first princi- 
ples for both AFM and FM ordering. In each case, 
as the frequency of the polar mode goes to zero at the 
corresponding rjc, the dielectric response diverges. For 
example, at 77=0.9%, we find 633^ 10'^ for the ground 
state AFM-PE phase. If a magnetic field is applied that 
aligns the spins into FM ordering, 633 increases 50-fold to 
that of the FM-PE phase, e33«5xl0''. Since the AFM- 
PE response remains finite as 77ci is approached while 
the FM-PE response diverges, the magnetocapacitance 
diverges in the vicinity of the phase boundary. 

Once the critical strain is crossed for each ordering, 
the unstable polar phonon freezes in and the ground state 
crystallographic symmetry becomes PAmm. In Fig.[2tb), 
we show the computed spontaneous polarization for the 
relaxed ground state structures as a function of epitaxial 
strain. A large F-rj coupling is evident for both orderings. 
For strains above 1.2%, the AFM state also become EE, 
but the range of strains over which the AFM-FE state 
is the ground state is extremely narrow, with a first- 
order transition to the FM-FE phase occurring at 77 « 
1.25%. At this phase boundary, the polarization in the 
FM phase is already over 10 fiC cm~^. Thus, we show 
that a previous unknown multiferroic phase of EuTiOs, 
with a symmetry-allowed linear magnetoelectric effect, 
can be stabilized by readily accessible epitaxial strains. 

The most interesting behavior in this system is pro- 
duced in the intermediate strain region, rjci < i] < 77^2, 



under applied electric £ or magnetic H. fields. We ob- 
tain a qualitative understanding of the phase diagram 
to leading order by adding the terms — M-Ti and —P £ 
to the energy, and taking M and P to have their zero- 
field values. For example, at 77=1.0% the ground state 
is still the AFM-PE phase but the application of a mag- 
netic field of sufficient strength to fully align the spins, 
©(IT) (the same order as that found by KT for the bulk), 
induces a substantial spontaneous polarization P«10/iC 
cm~^, by driving the system into the FM-FE phase. No- 
tably, P is several orders of magnitude greater than the 
polarization of all previously known EEs whose origin is 
through a coupling to spins. In addition, at 77=1.0% the 
application of an electric field, O{10^V cm"-'^) (small for 
thin films) favors the polar FE phase over that of the 
non-polar PE phase, inducing a magnetization of 7/i;, by 
driving the system into the FM-FE phase. We summarize 
the results for compressive strain in Fig. 13 

From these estimates, both magnetic control of the 
electric phase and electric control of the magnetic phase 
in epitaxially strained EuTi03 should be easily attain- 
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FIG. 2: (Top) Static dielectric constant £33 and (bottom) 
spontaneous polarization Pa as a function of epitaxial com- 
pressive strain rj. For £33, the solid lines are fits proportional 
to {ri-rjc)'^ while the dashed vertical lines indicate r/c. 
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FIG. 3: EuTiOa: Compressive epitaxial strain {r\) phase 
diagram. 



able experimentally. However, the symmetry of the PE- 
paramagnetic structure determines the lowest order spin- 
phonon coupling (the origin of the discussed effects) as 
a bi-quadratic effect Tint oc M^P^. So, while an ap- 
plied magnetic field can turn the polarization on from 
zero, there is no preferred direction and the magnetic 
field does not act to reversibly switch the polarization 
between symmetry- related orientations. At present, such 
switching has only been demonstrated in multiferroics 
with minute polarizations |lOl ll]| . Similarly, an applied 
electric field cannot be used to reversibly switch between 
different magnetization states [29|. though it can turn a 
substantial magnetization on from zero. 

In further analogy to SrTiOs, effects similar to those 
descibed above in EuTiOa are expected not only for 
compressive, but also for tensile epitaxial strain; first- 
principles investigations are in progress. The recent suc- 
cess in fabricating high-quality thin films of EuTiOa with 
zero epitaxial strain 30] is encouraging for experimental 
realization of the predicted effects. While we have fo- 
cused on the example of EuTiOa, we expect that addi- 
tional materials with the necessary characteristics could 
be identified and their epitaxial phase diagrams explored. 

In summary we have presented a realization of a new 
mechanism for strong coupling between magnetic and fer- 
roelectric ordering. The required material characteristics 
are a spin-phonon coupling through which FM spin align- 
ment softens a low-frequency polar mode that is strongly 
coupled to epitaxial strain. The predicted competition 
between an AFM-PE phase and a FM-FE phase allows 
magnetic phase control with an applied electric field, and 
electric phase control with an applied magnetic field, with 
modest critical fields. In addition to being a promising 
mechanism by which such phase-control can be achieved, 
we anticipate that epitaxial stabilization of a FM-FE 
ground state - such as that which occurs in EuTiOs above 
strains « 1.25% - over the bulk AFM-PE phase may also 
prove to be a useful avenue to the identification of new 
FM insulators suitable for spintronics applications. 
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